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We discuss electromagnetic modes in cavities formed by metamaterials with negative refraction and demon- 
strate that the straightforward approach to substitute negative values of the electric permittivity and the magnetic 
permeability leads to quasi-bound states with a negative quality factor To ensure positive quality factors and a 
consistent physical interpretation of the quasi-bound states it is essential to include the frequency dispersion of 
the permittivity and the permeability, as required by positive field energy and causality. The basic mode equa- 
tion and the boundary conditions including linear frequency dispersion are derived. As an example we consider 
a disk-like cavity with deformed cross sectional shape. The transition from the unphysical nondispersive case 
with negative quality factors to the dispersive case with positive quality factors is demonstrated numerically and 
in an analytical perturbative treatment. 

PACS numbers: 42.25.-p, 42.60.Da. 78.67.Pt 



I. INTRODUCTION 

Negative-index metamaterials (NIMs) are artificial com- 
posites characterized by simultaneously negative values of the 
electric permittivity e and the magnetic permeability /i |[l]-^ . 
These materials were theoretically predicted already in 1968 
by Veselago |4]. In present days there is a strong interest in 
such materials because of potential applications, such as elec- 
tromagnetic cloaking ||5[], sub wavelength imaging 16] and fo- 
cussing of light 

Cavities which confine electromagnetic waves in all three 
spatial dimensions have also attracted considerable attention 
in the recent years, in particular in the optical regime [SY The 
strong interest is partly due to the numerous future applica- 
tions, such as single-photon emitters and ultralow thresh- 
old lasers |[lOt4l2ll . and partly due to the possibility to ad- 
dress fundamental questions of light-matter interaction ifTsIl 
and ray-wave correspondence lfl4l - [l6ll . 

The fabrication of electromagnetic cavities made of NIMs 
is rather challenging with nowadays technology. The lit- 
erature in this field is therefore limited to a few theoreti- 
cal studies, e.g., on one-dimensional cavities made of dis- 
tributed Bragg reflectors 1 17] and two-dimensional superscat- 
terers [18]. These works, however, consider scattering of 
plane waves and do not treat the electromagnetic modes as 
quasi-bound states with finite lifetime. In this paper we fill 
this gap and discover a subtle difficulty when defining quasi- 
bounded states in the canonical way. We show that for a con- 
sistent physical interpretation of modes in NIMs the frequency 
dispersion of e and /i is crucial. Our consideration is gen- 
eral and applies to all sorts of NIM cavities. For illustration 
we present results for two-dimensional disk cavities with de- 



formed cross-sectional shape. 

The paper is organized as follows. Sections HIl and Ulllpro- 
vide a brief review on Maxwell's equations for monochro- 
matic fields and on electromagnetic cavities made of conven- 
tional materials. The boundary conditions for NIM cavites 
and the appearance of negative quality factors are discussed 
in Sec.|IV] Section [Vl deals with the basic properties of the 
frequency dispersion of e and /i. A modified mode equation 
including linear dispersion and a numerical solution of this 
equation is presented in Sec. [Vl] A discussion of the effects 
of the dispersion is given in Sec. I VIII Finally, Sec. lVIIll con- 
tains the conclusions. 



II. MAXWELL'S EQUATIONS 

The source-free Maxwell's equations in the frequency do- 
main are 

\/xE^i-B, VxH = -i-D, (1) 
c c 



V-L' = 0, V-B = 0, 



(2) 



where c is the speed of light in vacuum. The constitutive re- 
lations for monochromatic fields with frequency u! are in the 
isotropic and linear regime 



B{r,uj) fj,{f,uj)H{f,uj) 



(3) 
(4) 



with electric permittivity e and magnetic permeability /i. The 
boundary conditions at an interface between a material 1 and 
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a material 2 are given by 

Vy.{Ei-E2) - = i)y.{Hi- H2) , (5) 
iJ ■ {pL^H^ - PL2H2) = 0= (£1^1 -£2^2), (6) 

where v is the local normal vector on the boundary. 

III. CONVENTIONAL CAVITIES 

When discussing electromagnetic modes in cavities made 
of conventional materials one usually takes advantage of two 
simplifications: (i) the frequency dispersion of the permittiv- 
ity £ is ignored, assuming that the frequency interval of inter- 
est is sufficiently small, (ii) The permeability /i is assumed to 
be constant and unity at all frequencies, throughout the whole 
space. 

In the following, we focus on disk-like cavities mainly for 
illustration purposes. We emphasise here that our results can 
be extended to arbitrarily shaped three-dimensional cavities in 
a straightforward manner For the quasi-2D geometry of the 
disk one separates the (x, y)-dynamics in the plane of the disk 
from the z-dynamics by expanding the electric field in terms 
of E{x, y)e^^'^'^, with fc^ = or kz finite. In the latter case 
the refractive index n = y/ejl in the mode equation can be 
replaced by an effective index. Assuming a piecewise con- 
stant index of refraction one can derive the following mode 
equation for a quasi-two-dimensional disk 1, 19. .20,1 

-V^i^ = n^^^, (7) 

where fl is the frequency of the mode. This Helmholtz equa- 
tion holds for both transverse magnetic (TM) and transverse 
electric (TE) polarization. For TM polarization the elec- 
tric field is perpendicular to the cavity plane with = 
Re[^p{x, ?/)e~'^*]. At the boundary between a material 1 and 
a material 2 the general boundary conditions ^ and ^ give 
the following continuity relations for the wave function ip and 
its normal derivative di^tp along the normal i7 , 

^1 = ^2 , d,^i = d,iP2 (TM), (8) 

assuming that /ii = /i2. For TE polarization, the mag- 
netic field is perpendicular to the cavity plane with = 
Re[ip{x,y)e~'^^*]. The boundary conditions are 

'ipi ^'ip2 , -Kduipi = \d^'ip2 (TE), (9) 

again assuming jii — 112- At infinity, outgoing wave condi- 
tions in the two-dimensional disk plane 

with wave number k ~ 17 /c are imposed for both polariza- 
tions, which results in quasi-bound states with frequencies D, 
situated in the lower half of the complex plane. Whereas the 



real part is the usual frequency, the imaginary part is related to 
the lifetime r — — 1 / [2 Im il] . The quality factor of a quasi- 
bound state is defined by Q = —Re Vt/\21n\ Vl]. These reso- 
nant states, first introduced by Gamow [21] and by Kapur and 
Peierls L22,1 . are connected to the peak structure in scattering 
spectra; see ll23ll for an introduction. 

As an example we choose a disk-like cavity with the bound- 
ary curve being the limafon of Pascal, which reads in polar 
coordinates 

= + ecos0) . (11) 

For vanishing deformation parameter e this gives the circu- 
lar disk with radius R. We choose a deformed disk with 
e = 0.43. Exactly this geometry has been studied for con- 
ventional materials in the contextjjf directed light emission 
from microlasers, theoretically ll24tl as well as experimen- 
tally idiiim. The value of R itself is not relevant, only 
the ratio R/\ is important, where A 27r/fc is the wave- 
length. Therefore, we consider in the following a normalized 
frequency flR/c = kR. Figure [T] shows a typical TM polar- 
ized mode in a lima§on cavity with low index of refraction, 
n = 1.5, computed with the boundary element method ll28ll . 
The mode is localized along an unstable periodic ray trajec- 
tory, i.e., it is a so-called scarred mode 12911 . From the relative 
intensity of the different segments we can assess the direction 
of energy flow as indicated by the arrows. The frequency is 
riR/c = 44.9376 - i0.0622, the Q-factor is therefore about 
361. Such a medium-Q, scarred mode is well suited to demon- 
strate how the light is (partially) refracted out. This will be 
useful in the following comparison to the NIM cavities. 




FIG. 1: (color online). Near-field intensity pattern of an electromag- 
netic mode with normalized frequency Q,R/c — 44.9376 — i0.0622 
in a conventional cavity with e — 9/4 and fi = 1 (n — 1.5) sur- 
rounded by vacuum with e — /j. = 1. Arrows illustrate the direction 
of the energy flow. 



IV. NIM BOUNDARY CONDITIONS 

In the case of a NIM cavity the permeability ^ can no longer 
to be treated as spatially uniform as this quantity changes sign 
at the interface between of the NIM cavity and the surrounding 
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conventional material (in our case vacuum). In this situation 
the general boundary conditions Q and (|6]l give (see, e.g., 
Ref. Ill) 



V'l = ^-2 , —9.^1 = —d,iP2 (TM) (12) 

Ml M2 



instead of the special case in Eq. For TE polarization we 
have 



^1=^2, -dui^i = -a,V2 (TE) . (13) 
ei £2 



V. FREQUENCY DISPERSION 

When the dispersion of e and /i is important the expression 
for the field energy density ( fT4l ) has to be replaced by 



Stt V duj du! 



(15) 



The energy density W defined by Eq. (fTSb is positive provided 
that 



^>o 



did 



duj 



(16) 



Using these boundary conditions together with the 
outgoing-wave conditions (fTOl i at infinity we find that solu- 
tions of the mode equation (|7]i always have a negative Q-factor 
for NIMs, i.e., the intensity of such a solution ip does not de- 
cay but instead increases exponentially in time, which is un- 
physical for a passive material. As a typical example which 
will be discussed later in more detail we mention a mode in 
the limagon cavity with e = —9/4 and /i — —1 {\n\ — 1.5). 
The normalized frequency is VLqR/c — 45.1060 + i0.1647. 
The quality factor is therefore about Qq = —137. 

Note that for a conventional material a similar effect may 
occur if a thin active layer with very strong gain is placed 
at the interface. In that case the jump of the derivative im- 
posed for negative index materials in Eq. ( fTZt or (flJl l would 
be caused by an outflow of energy from that amplifying layer. 

It seems that for a passive NIM the exponential increase 
of the electromagnetic intensity contradicts the law of energy 
conservation. This is, however, not the case as can be seen by 
considering the electromagnetic field energy density 



Stt V 



(14) 



In this equation and also in the following ones we suppress 
the dependency on the spatial coordinates for notational con- 
venience. As already pointed out by Veselago |4], as e, /i < 
the field energy ( fT4l i in a NIM would be negative and un- 
bounded from below. In our case of a quasi-bound cavity 
mode this results in an exponential decay of the field energy 
towards — oo, as the cavity permanently looses (positive) en- 
ergy to the outside. Hence, the field intensity > Q inside 
the cavity increases exponentially in time. The imaginary part 
of the frequency 17 is therefore positive, and the quality factor 
negative, consistent with our numerical finding. 

It is worth mentioning that Dirac's wave equation for rel- 
ativistic electrons possesses a similar "radiation catastrophe" 
which disappears in a proper quantum field theoretical treat- 
ment. For the NIM materials, the problem of negative field 
energy density can be, however, solved already on the wave 
equation level, namely by the inclusion of the frequency dis- 
persion of the electric permittivity and the magnetic perme- 
abihty illMl. 



for any values of E'^ and H^. These inequalities imply lower 
bounds for the derivatives 



9£ 

du) 



> 



e 



9/i 



> 



UJ 



(17) 



with Lu > 0. For conventional materials with positive e and /i 
and nonresonant response to external fields the derivatives are 
small, and therefore often frequency dispersion can be safely 
ignored. For NIMs with negative e and fi this is never possi- 
ble. 

In a transparency region, where we can neglect absorption 
in the material, causality requires the following additional in- 
equalities 1I30I1 



de_ 2(1 - 

duj UJ 







^ > ^ 

did UJ 



m) 



(18) 



For negative e and fj, these inequalities are stronger than the 
ones in Eq.lfTTIi. Note that there is a controversy in the NIM 
community concerning the existence of negative refraction in 
such a transparency region; see, e.g., Refs. ll3ll - l33[l . This con- 
troversy is, however, not settled yet; so in the following we 
assume that a transparency region with negative e and ^ ex- 
ists. 

It is worth mentioning that in (unrealistic) materials hav- 
ing no absorption at all in the whole frequency interval, the 
inequalities in Eq. ( fTSl l can turn into equalities. This is, for 
instance, the case for the nonlossy Drude model with the di- 
electric function 



e{uj) = 1 



(19) 



and plasma frequency ujp. Another, more general model 
which is often used to describe NIMs locally in frequency 
space is the Drude-Lorentz system 



siuj) = 1 - 



170; 



(20) 



In the limiting case of no absorption 7 — > the inequali- 
ties in Eq. (fTSl l are fulfilled. For the subtle issue of causality 
in the nonlossy Drude-Lorentz model we refer the reader to 
Ref. la. 
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VI. MODIFIED MODE EQUATION 

In the following, we show that when dealing with quasi- 
bound states not only Maxwell's equations but also the con- 
stitutive relations including the frequency dispersion have to 
be extended to the complex frequency plane. This fact di- 
rectly leads to a modified mode equation which is capable of 
describing linear frequency dispersion. This mode equation 
predicts positive quality factors in agreement with the require- 
ments of positive field energy and causality. 

The key observation to start with is that a quasi-bound state 



(21) 



with complex-valued frequency is not a monochromatic 
wave as its Fourier decomposition 



En{t) = / dujE{u)e- 



(22) 



gives nonvanishing E{uj) in a region around oj ss Re(57). With 
this decomposition we can write the constitutive relation (|3]l 

as 



duj e{io)E{uj)e 



(23) 



To proceed further, let us fix a real-valued frequency Wr around 
which we study electromagnetic modes. We restrict ourselves 
to modes with sufficiently small |Im(51)|, i.e., not too small 
Q-factor, and Re (51) « cj^. This allows the expansion of the 
permittivity 



de 

e{uj) « e(wr) + t;- ~ ^i) 

OiU uir 



(24) 



Inserting this expansion into Eq 



gives 



e{LUr)Enit) + — 



^i—-ujr]Enit).(25) 



Now we exploit the property of quasi-bound states ^i?si 
—iilEn, which can be deduced from Eq. ( |2T]) and leads to 



Dnit) = s{n)Enit), 



with modified permittivity 



de 



(26) 



(27) 



Equation ( [26] l represents an analytic continuation of the per- 
mittivity e{uj) to the complex-frequency plane ft. For this 
continuation the linearization in Eq. (l24l i is not needed. In 
fact, for the modified mode equation that we derive in the fol- 
lowing, the extension to any analytic function e{uj) is straight- 
forward. Nevertheless, for clarity we restrict ourself to a lin- 
ear frequency dispersion. We would like to point out that the 
imaginary part of e is related to the frequency dispersion and 
not to optical absorption in the material, which is neglected 



here. Nevertheless, a realistic Q luy with negative imag- 
inary part introduces a kind of loss in the originally lossless 
medium which counteracts the nonphysical exponential in- 
crease of the intensity, thereby turning the negative quality 
factor into a positive one. 

In the following it will be convenient to express the values 
of the derivatives of e and /i at the fixed frequency Wr by their 
dimensionless hnear dispersions 



9£ 



e 



9/i 
duj 



(28) 



For NIMs (e < 0, /i < 0) these quantities a^, have to 
be chosen larger than 1 to satisfy the inequalities ( fTTI i (for 
conventional dielectrics a^, must be smaller than 1.). To 
satisfy the inequalities ( fTSl) for NIMs we must have 

2 2 

ae>2 - - , a^>2 . 

e fi 

With the quantities in Eq. ( |28] | we can rewrite Eq. ( |27] | as 



(29) 



e{n) = eiio,) 1 + 



UJr 



(30) 



We can derive a modified permeability in an analogue way 

ujy — Q 



(31) 



As a result of our considerations we can use Maxwell's 
equations ([Hi-© and the constitutive relations ©-(HI) for 
monochromatic waves with modified permittivity and perme- 
ability given by Eqs. ( [30l ) and ( [3T1 ). A direct consequence is 
the modified mode equation 



with the modified refractive index 



n{n) = y/e{n)il{Q) « n{ujr) 1 + 



and the dimensionless linear dispersion 



dn 



UJr 
UJr n 



(32) 



(33) 



(34) 



In the derivation we have ignored terms of order (Jl — Wr)^, 
which is consistent with Eq. ( l24b . For the square root in 
Eq. ( l33T l we choose the positive branch. Note that the (modi- 
fied) refractive index can be defined negative or positive. This 
does not matter for our purpose, as the sign of the refractive 
index neither enters the mode equation (l32b nor the modified 
boundary conditions 

V-i = ^2 , ^a.Vi = ^d,^2 (TM) (35) 

Ml 

V'l = ■02 , —dyipi = —dy%lj2 (TE) . (36) 
£1 £2 
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The phenomenon of negative refraction is here a result of the 
relative sign of the permittivity and the permeability in the 
boundary conditions ( l35l l and ( |36] |. 

After the quantities e, fi and their first derivative are spec- 
ified at a given reference frequency Wr, the modified mode 
equation (l32b . the modified permittivity (l30l l. the modified 
permeability dSTI ). the modified refractive index i33[ . and the 
modified boundary conditions ([35])-(l36l) have to be solved 
self-consistently. This can be done with only slight modifi- 
cations of standard approaches, such as the boundary element 
method iH. 

As an example we consider again a mode in the limagon 
cavity with normalized frequency Re(il)i?/c around 45. 1. We 
fix the reference frequency ui-^R/c therefore to be 45.1. Note 
that the precise value of the reference frequency is not rele- 
vant as long as we choose the values for the permittivity, the 
permeability and their derivatives accordingly, e.g. by reading 
off their values from a material dispersion curve at the given 
reference frequency. For the NIM we consider £(a;r) = —9/4 
and ii{LLir) — —1 as in the previous section. To fulfill the con- 
straints from Eq. ( |29] l 

g 

ae > 2 + - , > 4 , (37) 

we first set — — 4.1 and vary this value later. Accord- 
ing to Eq. ( [34l l the quantity q;„ is then also 4.1. Figure |2ja) 
shows as an example a TM polarized mode. The Q-factor 
turns out to be positive, Q ~ 438, due to the inclusion of the 
dispersion. Interestingly, the spatial mode pattern is rather in- 
sensitive to the inclusion of frequency dispersion, as a closer 
comparison of the dispersive case in Fig. |2ja) and the (un- 
physical) nondispersive case in Fig.|2|b) shows. 




FIG. 2: (color online), (a) Near-field intensity pattern of an electro- 
magnetic mode with ftR/c — 45.0966 — i0.0514 in a NIM cavity 
with e(ajr) = -9/4, /^(cjr) = -1 (\n{LOi)\ = 1.5), uJrR/c = 45.1, 
and linear frequency dispersion = — 4.1. Arrows illus- 
trate the direction of the energy flow, (b) Unphysical mode with 
QoR/c = 45.1060 + i0.1647 in a NIM cavity with e = -9/4 
and fi — —1 (|n| = 1.5) calculated without frequency dispersion. 

Contrasting the mode in the NIM cavity in Fig. 12 a) with 
a corresponding mode in a conventional material in Fig. [T] 
clearly reveals the negative refraction. To be more precise, 
one part of the beam is confined by total internal reflection 
and another part is refracted out. Moreover, a careful inspec- 
tion shows a different Goos-Hanchen shift (GHS) for the NIM 



and for the conventional material, as can be seen in Fig. [3] The 
GHS is a lateral shift of totally reflected beams along the op- 
tical interface due to interference jssll (for GHS in cavities 
see Refs. |36', 'jT]). According to Ref. pS], a light beam in a 
conventional material reflected at the interface to a NIM ex- 
periences a negative GHS. In our case the light propagates in 
the NIM and is reflected at the interface to a conventional ma- 
terial. We compute the GHS by reflecting a Gaussian beam 
at a planar dielectric interface, neglecting boundary curvature 
effects 1 39]. For the mode in the conventional cavity, see 
Fig. [3la), we find a nearly perfect agreement between such 
a beam reflection at a planar interface and the full mode cal- 
culation. For the mode in the NIM cavity, see Fig. Ob), the 
agreement is somewhat reduced, but nevertheless the appear- 
ance of a negative GHS can be clearly seen. 




FIG. 3: (color online). Goos-Hanchen shift in (a) a conventional cav- 
ity (cf. Fig.[TJ, and (b) a NIM cavity [cf. Fig.|2l;a)]. The dotted line is 
the center of a Gaussian beam being reflected at the dielectric inter- 
face. The appearance of the Goos-Hanchen shift along the boundary 
is evident. The mode pattern is the full numerical solution of the 
respective mode equation. 



VII. DISCUSSION 

To quantify small differences in the spatial pattern of the 
mode tpo without frequency dispersion [as, e.g., in Fig. |2|b)] 
and a mode i/; with dispersion [as, e.g., in Fig.|2|a)] we exam- 
ine the normalized spatial overlap 

S = \Jcd^dytP*tP\ ^^^^ 

^ dxdy ipo'ijjo'^ Jc dxdy 'ip*-)p 

We restrict the integrals to the interior of the cavity C as the 
exterior is influenced by the actual value of the quality factor. 
For the modes in Fig.|2](a) and (b) we find 1 — S* w 1.6 x 10~^, 
so indeed the overlap is nearly unity. The lower panel of Fig. |4] 
shows 1 — 5 as function of the linear frequency dispersion q;„. 
Except near = 1, the overlap 1 — S" is below 0.001. The 
upper panel of Fig. |4]shows Q/[—Qo) with Qq = —137 as 
function of «„. It can be observed that for Q!„ > 1, where the 
field energy is positive, the Q-factor is also positive. 

Can we understand the basic features observed in Fig. |4]? 
To do so, let us first note that the quantities |(n — n)/n\, 
|(e — e)/e|, and \ {fi — ^) / ^\ are small, as required by the linear 
expansions, e.g., in Eq. (|24] |. The applicability of perturbation 
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positive field energy ■ 



i-causality-5 




c/5 0.001; 



0.0001 r o 
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a = a 



FIG. 4: (color online). Q/{~Qo) (upper panel) and overlap dif- 
ference 1 — S [lower panel, cf. Eq.(l38b1 vs. linear frequency dis- 
persion a„ ~ cte = Q/j. All quantities are dimensionless. Here, 
Qo = —137 is the quality factor in the nondispersive case = 0. 
In the regime > 1 the electromagnetic field energy is positive 
and for > 4 also causality holds. The empty dots mark the re- 
sults of the full solution of the mode equation l |32| >. The solid line is 
the theoretical prediction in Eq. i39\ . 



theory to the Helmholtz equations (|7]i and ( l32b therefore im- 
pHes that the spatial mode pattern does not depend much on 
the frequency dispersion. But why is the Q-factor so strongly 
dependent on the frequency dispersion? This can be under- 
stood from the observation that for long-lived modes with high 
Q = — Reri/[2Imf2], the imaginary part of the frequency 
is small compared to the real part of fl. Therefore, even a 
small modification of Im(il) in absolute numbers can have a 
relatively large effect on the quality factor. 

To see how Q changes with the linear dispersion a„, con- 
sider the mode equation (Q with frequency fig and the modi- 
fied mode equation ( |32] | with frequency il. In both cases use 
the same negative e, /i. Both mode equations can give the 
same spatial mode pattern provided that nilo = nfi. From 
this relation we find for \oji- — Re(f2)| <C ojr and |Im(ri)| <C 
Re(ri) that 



Q_ 

-Qo 



(39) 



This shows that the modified refractive index turns the (non- 
physical) negative quality factor Qo into a positive one, Q > 
as soon as q;„ > 1, i.e., exactly under the condition which 
ensures a positive field energy in the NIM. To illustrate this 
relation from another point of view, let us rewrite the right 
hand side of Eq. ( |39] | as 



a„ - 1 = 



(40) 



where Vp — c/n is the phase velocity and Vg — dcu/dk is 
the group velocity with wave number k = nujjc. The in- 
equalities ( [TtT i required by a positive field energy ensure that 



the phase velocity and the group velocity Vg have a differ- 
ent sign in a NIM. The inequalities ( fTSl ) required by causality 
carry over to \vg\ < \vp\, i.e., superluminal energy propaga- 
tion is forbidden. 

As demonstrated in Fig. |4] the expression in Eq. ( [39l ) is 
in excellent agreement with the full solution of the modified 
mode equation ([32]). Note that near a„ = 1, where Q w 0, 
the absolute value of the imaginary part of and the quantities 
\{h — n)/n\, |(e — e)/e|, and — /i) are not small. Hence, 
the linear approximations in our theory are not justified, which 
explains why in this region 1 — 5 is larger. 

From another point of view, the expression in Eq. ( [39l ) can 
also be understood purely in real frequency space. Consider 
a resonant structure of some spectrum, let's say the Wigner 
delay time or a scattering cross section. Note that such res- 
onances are determined only by the value of the product 

= n{uj)oj of refractive index n and the frequency oj. In our 
example of the limagon cavity it is nkR rather than just kR, 
with wave number k = ut/c. Let the full width at half max- 
imum (FWHM) of the resonant peak be SC = nSuj + ujdn. 
Comparing the nondispersive case (Sn ~ 0) with the disper- 
sive case we get 



(41) 



again assuming that the change in the refractive index does 
not change the spatial mode structure. With Q = uj/Soj and 
Qo = oj/ SujQ we arrive after a few algebraic manipulations at 
Eq. 

A related expression as in Eq. ( [39] ) exists also for conven- 
tional materials where it was used to predict enhancement of 
quality factors in microcavities using highly dispersive mate- 
rials |40]. For small group velocity Vg, corresponding to slow 
light, the effect of the Q-factor enhancement in Eqs. (39[ and 
( [40] ) is strongest. This prediction has been confirmed recently 
in experiments on slow light in photonic crystals |41]. 

Finally, we note that scattering of a monochromatic wave 
with (real-valued) frequency uj at an obstacle made of a NIM 
is described by the ordinary mode equation ([TJ using incoming 
and outgoing wave conditions |17, 18]. However, when cal- 
culating spectra it is necessary to use the frequency-dependent 
permittivity e{uj) and permeability to be consistent with 
the requirement of causality. This is usually ignored, perhaps 
because in a spectrum the FWHM Suj cannot be easily distin- 
guished from —duj. 



VIII. CONCLUSIONS 

We addressed quasi-bound electromagnetic modes in 
negative-index metamaterial cavities. The simple approach 
which substitutes negative values of the electric permittiv- 
ity and magnetic permeability into the boundary conditions 
(which suffices to obtain negative refraction) gives rise to 
modes with negative Q-factor. This unphysical behavior can 
be removed by including linear frequency dispersion in the 
mode equation and in the boundary conditions, as is required 
by a positive field energy and causality. At complex reso- 
nance frequency, the effective permittivity then acquires a fi- 
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nite imaginary part, which attenuates the mode even in ab- 
sence of physical absorption. 

As an example we studied a disk-like cavity with noncircu- 
lar cross-sectional shape. The modified mode equation results 
in modes with positive quaUty factor and clear signatures of 
negative refraction and negative Goos-Hanchen shift. 
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